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A method to provide active noise control to reduce noise and 
vibration in reverberant acoustic enclosures such as aircraft, 
vehicles, appliances, instruments, industrial equipment and 
the like is presented. A continuous-time multi-input multi- 
output (MIMO) state space mathematical model of the plant 
is obtained via analytical modeling and system identification. 
Compensation is designed to render the mathematical model 
passive in the sense of mathematical system theory. The com- 
pensated system is checked to ensure robustness of the pas- 
sive property of the plant. The check ensures that the passivity 
is preserved if the mathematical model parameters are per- 
turbed from nominal values. A passivity -based controller is 
designed and verified using numerical simulations and then 
tested. The controller is designed so that the resulting closed- 
loop response shows the desired noise reduction. 
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METHOD AND SYSTEM TO PERFORM 
ENERGY-EXTRACTION BASED ACTIVE 
NOISE CONTROL 

STATEMENT REGARDING FEDERALLY 5 

SPONSORED RESEARCH AND DEVELOPMENT 

This invention was made in part with Government support 
under Grant Number NCC-1-01039 awarded by NASA and 
Grant Numbers 0196198 and 0301740 awarded by the N SF. l o 
The Government may have certain rights in this invention. 

FIELD OF THE INVENTION 

This invention pertains to noise suppression, and more 15 
particularly, relates to active noise control. 

BACKGROUND OF THE INVENTION 

Sound is a longitudinal wave phenomenon. Interference 20 
between waves has been known and observed since the time 
of Sir Isaac Newton. Active suppression of noise has been a 
goal for many decades. The potential applications of active 
noise suppression are many, ranging from aircraft cabins to 
automobile interiors to medical equipment (e.g., dentist 25 
drills, operating rooms, etc.) to buildings to household 
machinery such as vacuum cleaners, lawn mowers, snow 
blowers, washers, dryers, etc. The concept of wave dynamics 
and the ability of electronic amplifiers and moving coil speak- 
ers to generate accurate sound waves led to techniques used in 30 
the 1930s to reduce noise for a harmonic acoustic distur- 
bance. The next important step was taken in the 1950s with 
the experimental demonstration of a very high-gain static 
feedback acoustic noise reduction scheme. In today’s termi- 
nology, it is called a collocated control scheme with its inher- 35 
ent robustness properties. This basic collocated high-gain 
static feedback scheme is now commercially used for indus- 
trial ear protection devices. 

In the 1970s, the development of acoustic noise control 
took a slightly different turn. Papers such as those by Jessel & 40 
Mangiante, and Swinbanks addressed the problem of acous- 
tic noise control in ducts. The target application discussed in 
these papers was the reduction of noise in HVAC systems. 
The main idea was to measure the noise disturbance using a 
microphone located at an upstream location and then to feed- 45 
forward this disturbance via a controller to a speaker located 
somewhere at the downstream location. The research showed 
that a proper adjustment of the controller gains would arrest 
the propagation of the disturbance downstream from the con- 
trolling speaker. This control technique, referred to as the 50 
feedforward control, led to several years of active noise can- 
cellation research and development using the feedforward 
approach. 

For several years, feedforward methods were very popular. 

It was argued that since feedforward scheme doesn’t alter 55 
system dynamics, there was no risk of destabilizing the sys- 
tem (or making it any worse). These schemes required (or 
otherwise assumed) that the disturbance measuring micro- 
phone (located upstream) does not sense the output from the 
control speaker (located downstream). If it did, then it would 60 
be a feedback scenario with a danger of closed-loop instabili- 
ties. This meant that the schemes could be successfully imple- 
mented only with unidirectional microphones and speakers. 

In general, these conditions are very stringent. Nevertheless, 
researchers came up with ingenious arrangements and 65 
obtained impressive results. Most well-known are the experi- 
ments based on feedforward control by Ross and Roure 
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wherein the controller was designed to perfectly cancel the 
acoustic noise. These experiments involve exact knowledge 
of several system transfer functions, i.e., the design needs to 
calculate differences of transfer lunctions and invert them 
leading to high sensitivity of the controllers to uncertainty in 
the system model. Owing to high sensitivity, the model infor- 
mation has to be obtained experimentally. This remedy, how- 
ever, did not alleviate the problems as the controllers were 
very sensitive to the implementation and finite digit arithmet- 
ics and had to be implemented with on-line adaptation 
schemes. In spite of these hurdles, their success was a big 
inspiration for much of the work in active control of acoustic 
noise. 

However, one important drawback of the prior feedforward 
strategies that are used is the necessity of an on-line adapta- 
tion, not due to the nature of the system (acoustic duct) but 
due to the nature of the controller itself. Adaptive control 
comes at a price. For large bandwidth systems such as a 3-D 
acoustic duct, support of a fast digital signal processor with 
significant memory is essential for the feedforward tech- 
niques to work. Therefore, the implementation of feedfor- 
ward controllers really got boosted only after the arrival of 
cheap, high performance DSP (digital signal processor) 
chips. Another difficulty with feedforward schemes is to 
ensure that the adaptive control algorithm will always con- 
verge. Most algorithms guarantee convergence with known 
model, known controller order, and known persistency of 
excitation but, unfortunately, such ideal conditions are very 
hard to meet. As a result, feedforward methods for noise 
control are generally effective only at specific locations, usu- 
ally downstream of the disturbance, and at specific distur- 
bance frequencies. The need to use high performance DSPs 
for the computationally intensive on-line adaption schemes 
results in an expensive solution and the lack of robustness can 
make the noise worse in the event of small deviations. 

Since the late nineteen eighties, active noise control 
research took a multifaceted approach. During this time, 
research efforts extended to a multitude of problems, includ- 
ing such items as investigation of structure-borne noise, 
acoustic radiation from structures, reduction of noise in ducts 
and enclosures, and acoustic -structure interaction. Feedfor- 
ward adaptive algorithms work well when a signal strongly 
correlated with the noise source is available. Most of the 
algorithms converge nicely if the dynamics of the path 
between the control actuator and the sensor is known. 
Although some earlier successful solutions to structure-borne 
sound used multiple microphones and loudspeakers, subse- 
quent studies showed that loudspeakers alone are not very 
efficient in controlling structure-borne sound. This observa- 
tion led researchers to invent new actuators that directly act on 
the structure itself. It was found that an effective way to 
reduce the radiated acoustic energy from the structure is to 
control the structure so as to minimize the radiated acoustic 
power. Although a majority of the research during the 1 980’ s 
and 1 990’s still used feedforward adaptive algorithms, a large 
number of researchers were also engaged in developing and 
using feedback techniques. 

The research in the nineteen nineties especially involved 
fundamentally new approaches to active noise control — a 
feedback control instead of traditional feedforward methods. 
A significant effort was also expended looking at specific 
problems (e.g., active noise control for aircraft cabins, control 
of structure-borne noise, and a few others). In spite of diffi- 
culties, a number of approaches have succeeded in feedback 
control of reverberant enclosures but they are restricted pri- 
marily to the low-frequency bands. Additionally, much of the 
experimental work was restricted to small size cavities. Nev- 
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ertheless, the success stories of feedback methods in nineties 
gave boost to the researchers for further exploration of feed- 
back concept for active noise control. 

The major hurdles in successful implementation of feed- 
forward methods included availability of a reference signal 5 
having strong correlation to the noise source, necessity of 
online adaptation, and lack of controller robustness. For 
example, for structure-borne sound in aircraft it is difficult to 
obtain a reference signal which is strongly correlated with the 
noise producing mechanism. This makes feedforward control to 
very ineffective. The advances in robust control methodolo- 
gies coupled with knowledge gained in feedforward control 
schemes made it possible to develop closed-loop stable feed- 
back control methods. Moreover, feedback methods have 
natural robustness compared to feedforward methods which 15 
is inherent in the architecture. A comparison of LQG (linear- 
quadratic-Gaussian), rate feedback, and filtered-x LMS (least 
mean squares) algorithm, presented for SISO (single output, 
single input) feedback control of structure-acoustic dynam- 
ics, showed that there is a definite rationale to work with 20 
feedback methods. 

However, designing controllers for acoustic systems is 
very challenging since there is no natural roll-off at high 
frequencies and the systems are modally very rich. In the mid 
to late nineteen nineties, noteworthy research demonstrating 25 
effectiveness of feedback control methodologies emerged. 

For example, a PDE-based (partial differential equation- 
based) feedback control design was employed for controlling 
acoustic noise in 3-D enclosure using piezo-actuated struc- 
tures . The results were only numerical and the noise was only 30 
tonal. 

Although the feedback control methodology has great 
potential, it has several obstacles to overcome. One problem 
with many feedback control techniques used is that the noise 
attenuation can be guaranteed only locally (i.e., only at the 35 
sensor locations) and not uniformly in space and/or fre- 
quency. In fact, the noise levels can go high at some other 
locations although there is attenuation at the sensor locations. 
This essentially is the result of redistribution of acoustic 
energy in space. It is also possible to have worsening of noise 40 
levels at the same location due to redistribution of acoustic 
energy in frequency, commonly referred to as the waterbed 
effect. 

BRIEF SUMMARY OF THE INVENTION 45 

The invention provides a method and apparatus to provide 
active noise control to reduce noise and vibration in reverber- 
ant acoustic enclosures such as aircraft, vehicles, household 
appliances, high precision instruments requiring vibration/ 50 
noise isolation, industrial equipment and the like. The inven- 
tion allows the designs of quieter aircraft fuselages, vehicles, 
vacuum cleaners, washers, dryers, lawn and garden equip- 
ment (e.g., power blowers, lawn mowers, snowblowers, etc.), 
and other industrial/household machinery such as grinding 55 
machines, milling machines, cooling systems, and the like. 
The apparatus uses compensators and feedback control to 
render the system “passive” and reduce the structural and 
acoustic energy in the enclosure. 

The method to render the apparatus includes the steps 60 
described below. A continuous-time multi-input multi-output 
(MIMO) state space mathematical model of the plant (i.e., the 
acoustic enclosure) is obtained. The model is derived via 
analytical modeling and system identification. Compensation 
is designed to render the mathematical model passive in the 65 
sense of mathematical system theory. The compensation may 
be series, feed-forward, feedback, hybrid, etc. 
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The compensated system is checked to ensure robustness 
of the passivity property of the plant. The check ensures that 
the passivity is preserved if the mathematical model param- 
eters are perturbed from nominal values. A passivity -based 
controller is designed and verified using numerical simula- 
tions and then tested. 

Other aspects and advantages of the invention will become 
more apparent from the following detailed description when 
taken in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. lisa block diagram generally illustrating an acoustic 
environment in which the present invention operates; 

FIG. 2 is a flow chart illustrating the process of creating a 
controller design in accordance with the teachings of the 
present invention; 

FIG. 3 is a block diagram of an acoustic enclosure on which 
the controller design of the present invention operates; 

FIGS. 4a-4e illustrate passification methods according to 
the teachings of the present invention; 

FIG. 5 illustrates a passivity based controller structure in 
accordance with the present invention for a one-dimensional 
enclosure system; 

FIG. 6 is a graph illustrating experimentally obtained fre- 
quency-domain noise suppression response for the enclosure 
of FIG. 5 for a broadband disturbance; 

FIG. 7 is a graph illustrating noise suppression using the 
teachings of the present invention in a three-dimensional 
reverberant structure for a three-tone disturbance noise at 
frequencies of 100 Hz, 2256 Hz, and 290 Hz; 

FIG. 8 is a graph illustrating experimentally obtained noise 
suppression response for the system of FIG. 7; and 

FIG. 9 is an isometric view of a three dimensional enclo- 
sure having a vibrating surface. 

DETAILED DESCRIPTION OF THE INVENTION 

The invention extracts acoustic energy and structural 
energy in acoustic enclosures. The system utilizes acoustic 
sensors and actuators (e.g., microphones and loudspeakers) 
and actuators (e.g., accelerometers and piezoelectric devices) 
to reduce the acoustic energy and structural energy of the 
system. A compensator is designed to render the open-loop 
system passive in the sense of mathematical system theory. A 
positive-real feedback controller is then designed. The result- 
ing system is typically a multi-input multi-output feedback 
control system that incorporates the passifier and the positive- 
real controller. The system can achieve uniform (in space and 
frequency) broadband reduction in noise and vibration levels 
in the enclosures. As a result, quieter aircraft fuselages, 
vehicles, vacuum cleaners, washers, dryers, lawn and garden 
equipment (e.g., power blowers, lawn mowers, snow blowers, 
etc.), and other industrial/household machinery such as 
grinding machines, milling machines, cooling systems, and 
the like can be designed. 

Turning to the drawings, wherein like reference numerals 
refer to like elements, the invention is illustrated as being 
implemented in a suitable environment. FIG. 1 illustrates the 
control system 100 and acoustic enclosure 102. The acoustic 
enclosure 102 may be any of the acoustic enclosures 
described above (e.g., enclosures such as aircraft, vehicles, 
household appliances, high precision instruments requiring 
vibration/noise isolation, industrial equipment and the like. 
The control system 100 includes acoustic pressure sensors 
104 (e.g., microphones) and structural vibration sensors 106 
(e.g., accelerometers. The acoustic pressure sensors 104 are 



US 7,623,993 B2 


5 

placed within the acoustic enclosure 104 and the structural 
vibration sensors 106 are placed on the walls of the acoustic 
structure 104. Any number of acoustic pressure sensors and 
structural vibration sensors may be used. The location of the 
sensors 104, 106 may be application dependent. 

The sensors 104, 106 are in communication with signal 
conditioners 108, which consist primarily of analog filters. 
The conditioned outputs of the signal conditioners 108 are 
converted to discrete-time digitized signals by analog-to- 
digital (A/D) converters 110 located within energy -extraction 
controller 112. The energy-extracting controller 112 consists 
of a processing unit that stores data and computes the control 
(or command) signals. In one embodiment, the processing 
unit is a digital signal processor. The energy -extracting con- 
troller also includes a variety of computer readable media. 
Computer readable media can be any available media that can 
be accessed by controller 112 and includes both volatile and 
nonvolatile media, removable and non-removable media. By 
way of example, and not limitation, computer readable media 
may comprise computer storage media and communication 
media. Computer storage media includes both volatile and 
nonvolatile, removable and non-removable media imple- 
mented in any method or technology for storage of informa- 
tion such as computer readable instructions, data structures, 
program modules or other data. Computer storage media 
includes, but is not limited to, RAM, ROM, EEPROM, flash 
memory or other memory technology, CD-ROM, digital ver- 
satile disks (DVD) or other optical disk storage, magnetic 
cassettes, magnetic tape, magnetic disk storage or other mag- 
netic storage devices, or any other medium which can be used 
to store the desired information and which can be accessed by 
controller 112. Communication media typically embodies 
computer readable instructions, data structures, program 
modules or other data in a modulated data signal such as a 
carrier wave or other transport mechanism and includes any 
information delivery media. By way of example, and not 
limitation, communication media includes wired media such 
as a wired network or direct-wired connection, and wireless 
media such as acoustic, RF, infrared and other wireless 
media. Combinations of the any of the above should also be 
included within the scope of computer readable media. 

Digital-to -analog (D/A) converters 114 are embedded in 
the controller 112 and convert the digital control commands 
to analog commands. The analog commands are amplified by 
amplifiers 116. The amplified commands are sent to actuators 
118 to extract energy from the acoustic enclosure 102. The 
actuators 118 consist of one or more acoustic signal genera- 
tors 120 (e.g., loudspeakers) at various locations within the 
acoustic enclosure 102 and/or one or more piezo actuators 
122 attached to the inside or outside of the enclosure wall. 
While FIG. 1 shows the amplifiers connected to only one 
acoustic signal generator 120 and one piezo actuator 122, it is 
recognized that the amplifier is in communication with all 
active acoustic signal generators 120 and piezo actuators 122. 
The actuators 118 generate the commanded structural and/or 
acoustic wave motion. Performance monitoring sensors (not 
shown) may be part of the system to measure noise (i.e., 
pressure) and vibration (e.g., acceleration) at desired loca- 
tions to monitor the performance of the control system. The 
performance sensors consist of one or more acoustic pressure 
sensors (e.g., microphones) placed within the acoustic enclo- 
sure, and/or one or more structural vibration sensors (e.g, 
accelerometers) placed on the walls of the structure. 

Turning now to FIG. 2, the controller design extracts the 
acoustic energy and the structural energy in the system. The 
controller design is based on passivity theory. The controller 
is designed according to the following steps. A continuous- 
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time multi-input multi-output (MIMO) state-space math- 
ematical model of the acoustic enclosure (henceforth also 
called the “plant”) is obtained (step 200). The model is 
obtained via analytical modeling and system identification. 
5 Compensation to render the mathematical model “passive” in 
the sense of mathematical system theory is designed (step 
202). The compensation may be series, feed-forward, feed- 
back, hybrid, etc. (as described herein below). 

The passivity of the compensated system is checked (step 
10 204) using experimentally obtained frequency response data . 
Numerical tests are performed to check robustness of the 
passivity property of the plant, i.e., to ensure that the passivity 
is preserved if the mathematical model parameters are per- 
turbed from their nominal values. If this is not the case, 
15 compensation is redesigned (step 206). Steps 204-206 are 
repeated until the passivity is preserved. A passivity-based 
controller is designed using techniques described herein (step 
208). Numerical simulations are performed of the closed- 
loop system in the presence of a simulated broadband distur- 
20 bance input (step 210). If the performance is not satisfactory, 
the controller is redesigned and steps 208-210 are repeated 
until performance is satisfactory. 

Turning back to step 200, it should be noted that many 
members of industry who are engaged in active noise control 
25 research resort to experimental identification of system mod- 
els based on frequency response data. A typical procedure 
involves first obtaining frequency response data and then 
using this data to obtain necessary transfer function or state- 
space models. These identification techniques also have limi- 
30 tations as the currently available system identification algo- 
rithms work well only for low-order systems and systems 
with nicely separated modes. For a 3-D acoustic enclosure 
with a vibrating boundary surface, these identification tech- 
niques face difficulty even with two-input one-output system. 
35 The invention provides a general purpo se analytical modeling 
methodology for a 3-D acoustic enclosure having a vibrating 
boundary surface. This modeling technique accounts for the 
acoustic -structure interaction dynamics and yields time-do- 
main as well as frequency -domain models. 

40 A non-homogeneous wave equation is derived below for a 
general 3-D enclosure wherein the source term in the equa- 
tion represents the acoustic disturbance arising from the 
vibrating plate boundary. Since a closed-form analytical solu- 
tion is not possible when the entire boundary is considered to 
45 be vibrating, an assumption is made where only part of the 
boundary is considered to be vibrating. First, an infinite 
dimensional model is obtained which accounts for the struc- 
tural-acoustic coupling. Next, the assumed modes method is 
employed to obtain a finite-dimensional approximate model, 
50 which is used for control design purposes. Both the state- 
space model and the transfer function model are derived 
analytically. 

Turning now to FIG. 3, a 3-D acoustic enclosure 300 with 
a part 302 of one boundary surface 304 representing a vibrat- 
55 ing plate (denoted by S x ). Let the entire boundary surface be 
denoted by S. Prior to proceeding, given below is the list of 
definitions of various variables and parameters used in deriv- 
ing the acoustic -structure interaction dynamics for the enclo- 
sure. 

60 

v (t, x, y, z) is fluid velocity in m/s 

p(t, x, y, z) is acoustic pressure fluctuation from equilibrium 
in N/m 2 

p is fluid density fluctuation from equilibrium in kg/m 3 
65 p 0 is acoustic pressure at equilibrium in N/m 2 
p 0 is fluid density at equilibrium, kg/m 3 
co(t, x l5 y x ) is the transverse displacement of the plate in m 
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f(t, x l5 y x ) is the externally applied force on the plate due to 
acoustic pressure from an exterior noise field in N 
u z (t) are control forces acting at the points (x /f ; y /z ); i=l, 

2, . . . , r on the plate and normal to the plate in N 
The acoustic pressure and the fluid density at any instant of 5 
time are given by p(t, x, y, z)+p 0 and p(t, x, y, z)+p 0 . The 
control forces, u z (t), are produced by piezoelectric actuators 
and can be treated as point sources. For theoretical derivation, 
we make a necessary and reasonable assumption that the plate 
is simply supported on its four sides and that the initial dis- to 
placements and velocities at all points on the plate are zero. 
The goal is to obtain a governing relation between the acous- 
tic pressure fluctuation p(t, x, y, z) inside of a 3-D enclosure 
such as enclosure 300 in FIG. 3 and external disturbances f(t, 
x 1? y x ) and u z (t), i=l, 2, . . . , r. The model of the vibrating 15 
boundary surface (S x ) is first established followed by the 
development of the acoustic model for inside of the enclosure. 

The dynamic model of the vibrating surface 304 (i.e., the 
plate) can be obtained using Hamilton’s principle, which 
states that: Of all possible paths of motion that can be taken 20 
between the instants of time t l and t 2 , the actual path taken by 
the system gives a stationary (extremum) value to the integral 
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where h, and p^ are the thickness and density of the plate, 
respectively, E is the Young’s modulus, v is the Poisson’s 
ratio, and 
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is the flexural rigidity of the plate. In Equation 2, the function 
F is given by 
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where, K l is the kinetic energy of the system, U 1 the potential 
energy of the system, and W x is the work of the externally 30 
applied forces. It is assumed that the quantity K 1 +W 1 -U 1 is 
conservative. Now, for a two dimensional system with one 
dependent variable, if the integral in Equation 1 can be 
expressed as 



Ui(t)(o(t, xi.yMxi - xu) 


<50>i - yu) + [/(*> xu yO + p(t, xu yu 0)Mr, x u yi) 


The governing equation of the plate can be obtained from 
Equation 3 as 


/(«) = 


(2) 


1 3^1 5 7, OJ, OJ, OJ x l, OJyl , , OJ x ^y^ , OJy^y^ ) 

i 

dx i dyidt 


40 


4 d 2 o) 

DV oj + hp p -jj- = 


(6) 


^ J u i (t)6{x l -x li )S(y l -yu) + f(t,xi,yi) + p(t, xu yu 0) 


from Hamilton’ s Principle and variational calculus, the gov- 45 
eming PDE (partial differential equation) for the system can 
be derived from 


h (oj) ~h{6j) = 0 
where 


(3) 


k(oj) = F u - ■ 


dF t 


dxi 


dF^, 


to y i 


a 2 F,„, 


^ x l x l 


d 2 F [ ’ , 


dyi dx\ 

i /■ \ dF * 

h(l0) = -aT 


dyl 


d 2 Fu 

dxidyi 


50 
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The boundary conditions for the plate are given by 


d 2 oj(t, 0, y) 



dx 2 

oj(t, 0, _y) = 0 

d 2 oj(t, l x , y) 

oj(t, l x ,y) = 0 

dx 2 

oj(t, x, 0) = 0 

d 2 oj(t, x, 0) 

Oj(t, x,l y ) = Q 

dy 2 


d 2 Oj(t, X, ly) 


dy 2 


= 0 
= 0 
= 0 
= 0 


(7) 


For a vibrating plate, one can calculate U l5 K 1; and W x as 
follows: 


Vi = 


(4) 


D 

2 



d 2 oj 
dx 2 



-2(1 - v) 


d 2 oj d 2 oj 
dx\ dy\ 


( d 2 oj ) 2 1) 

Mr* 


The governing equation for the acoustics in the enclosure 
400 shall now be derived. For the enclosure 300 under con- 
sideration, a part of one of its boundary surfaces is a vibrating 
60 surface and we can consider this surface as a distributed noise 
source. The fluctuating volume flow per unit volume is then 
given by 

p 0 oo/r, x v _y 1 )6(z) for O^x^b-a and 0^y^d-c 

65 The linearized mass conservation equation and momentum 

conservation equation for the system under consideration can 
be written as 
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dp{t, x , y, z) 


+ p 0 V ■ v(t, x, y, z) = Po^ f (r, x u >i)<5(0 


d v (r, x, v, z) 

Po ^ + Vp(r, x, y,z) = 0 


Similarly, the equation to determine Y is 

} < 4 )_2(5 2 r-[(^:Q) 4 -p 4 /7=0 (17) 

For the case, (KQ) 4 -(3 4 >0, the solutions of Equations 16 
5 and 17 are given by 

X(x 1 )=A 1 cos(px 1 )+^ 2 sind^i) (18) 


Now, using the relationship p(t, x, y, z)=c 0 2 p(t, x, y, z), 
where c 0 is the sound speed, the divergence (V*) of the 10 
momentum conservation equation, and then subtracting the 
time derivative of the mass conservation equation, the non- 
homogeneous wave equation is obtained as follows 


7(y 1 )=5 1 cos(8y 1 )+5 2 sin(8y 1 )+5 3 cos l^ejY+i^ sin 
h(^i) 


<5 = ^(KfL) 2 -^ , s = y (Kft) 2 + f£ , 


V 2 p(t, x , y, Z)~ 

1 d 2 p(r,jt, y,z) =-po<*>tt(t,x-a,y-c) 

4 5? 

[H{x-a)-H{x-b)\ 
[H(y - c) - H(y - d)]d(z) 


and A z , B z are coefficients to be determined. Note that for the 
case when (KQ) 4 -(3 4 <0, there are no solutions. The solution 
W(x 1 ,y 1 ) can now be given by 


W(*i, yi) = [Aicos(/5*!) + A 2 sin08*i)] x 


where H( # ) is the Heaviside function. For the closed 3-D 25 
enclosure 300 under consideration we have the following 
boundary conditions 

v (t, x, y, z)' n =0 on S (1 1) 


[i?icos(<5^i) + Z?2sin(<5;yi) + 
^coshfe^i) + B4sinh(£_yi)] 


The boundary conditions (7) can be expressed in terms of 
Was 


where n is the outward normal direction at every point on 
surface S. Using Equation 9, the boundary conditions in Eq. 
(11) can be rewritten as 

Vp(t, x, y, z )• n =0 on .S' (12) 


r W(0, Vi) = 0, 
W(b- a, >i) = 0, 
W(xi, 0) = 0, 
,W(x l ,d-c) = 0, 


d 2 W( 0, yi ) 


d 2 W(b - a , yi) 


We will now derive the governing differential equations for | W(xi , 0) = 0, d 2 W(x { , 0) 

the vibrating surface (plate) as well as the acoustic pressure ( W(xi, d-c) = 0, dyj - 0 

fluctuation dynamics in the modal space as it will allow us to ^ d-c) 

obtain the finite-dimensional approximate models of the sys- 40 ~ = 0 

tern for controller design and experimental verification. First, 
we will consider only the homogeneous part of Equation 6 

with boundary conditions (i.e., Equation 7), that is, The coefficients and parameters in Equation 20 can be 

determined from these boundary conditions to be 
45 A 1 =B 1 =B 3 =B 4 =0, sin[(b-a)p]=0, sin[(d-c)8]=0. The normal 
nv 4 -i- h - n ^ modes of the plate are 


DV M + h Pp —= 0 


<b nm (x i, yi) =A nm sin(p rt x 1 )sin(5 M r 1 ), n, m=l, 2, 


Assume the solution is of the form co(t, x l5 y 1 )=Re{W(x 1 , 50 where, (3 w =n:t/(b-a), 6 m =rnjt/(d-c). The natural frequencies 
yi)e ia *}. Equation 13 is transformed into are 


V 4 W- (KQ,) 4 W=0 


where K^hp^/D, Q 4 =co 2 . Now separating variables x : and 
y l5 W can be written as 

Wix^Xix^Yiy o 

Equation 14 can now be re-written as 

[Y 4) - (KafX]Y+2X"Y"+XY ( *>=Q (1 5) 

Now considering the solutions which are separable. Let 

X"=-$ 2 X (16) 


' K 2 i\b- 


n \2 f m \2n 


Y 4) =-pTx"=p 4 x 


By the orthogonality condition: 


•S',/ 


<pn 2 m 2 dS ! = 6{n u n 2 )S(m l , m 2 ) 


65 the coefficients can be determined as A nm =2. 

The excitation of the plate due to acoustic pressure load 
from inside of enclosure is insignificant compared to piezo 
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11 

forces. That is, the dynamic interaction is much stronger from 
structure to acoustics than from acoustics to structure. Data 
obtained in applications such as aircraft cabins shows that the 
dynamic interaction from fuselage to interior acoustics is 
much stronger than in the opposite direction. In view of this, 5 
the acoustic pressure fluctuation p(t, x, y, z) term in equation 
6 is neglected without much loss of accuracy. The mathemati- 
cal simplicity achieved by this assumption outweighs the 
effects of modeling errors caused by the omission of this term. 

As a result, the governing equation 6 for the plate can now be 
re-written as 


. V CO V— 1 

DV 0J + h PfJ^T = Aj u ' 


= u i (t)6(x l -xu] 


(24) 


15 


f(t,x i, yO 


Using the assumed modes approach the solution of this 20 
PDE can be expressed as an infinite series representing an 
infinite sum of the weighted normal modes of the system, i.e., 


U)(t, *1 , >i) = ^ ^ qrmiWnmix i, }>i) 
n=l m= 1 


where q wm (t) are the modal coordinates and (|) wm (x 1 , y : ) are the 30 
corresponding mode shapes. For the control design purposes, 
we need to obtain a finite dimensional model of the system. 
This can be done by approximating the solution in equation 
25 using a finite number of modes. That is, we can truncate the 
series by choosing limits on n and m as n: 1 — >p : and m: 1 — >p 2 . 35 
Then, the solution for co(*) can be given by 


pi pi (26) 

0 )(t, Xl , > ! ) « ^ ^ (*1 , ?! ) 40 


Substituting equation 26 into equation 24 and multiplying 
it by ^(x^ yj and integrating it over S l5 we get 


bppQnmW + D (P 2 n + 8 lf £«m( 0 = ^ Y. «/ (O0 nwl (*i;, Jli) + 

fnmit) 

where, 

frm(t) = /h, Xl, yi)<pnm(Xi,yi)dSL 

It should be noted that equation 27 does not contain any 
damping terms as there is no known method to analytically 
derive the modal damping in the system. The damping has to 
be introduced a posteriori into the model. Introducing the 
modal damping £ wm (n=l, 2, . . . , p 1; m=l, 2, . . . , p 2 ) into the 0 
equation, we obtain 


45 

(27) 

50 

55 


q nm (t) + '2^ nm o) nm q mn {t)+ i 

ulmqnmit) h Pp 


(28) 


-continued 



«iWim(*l. }h) + frmAt) 


The modal coordinates can be determined by Duhamal inte- 
gral. Once the modal coordinate solution is obtained, the 
response of the plate to external and control forces can be 
obtained by equation 26. 

The solution for acoustic pressure p(*) using the assumed 
modes method. Let us assume that the solution of the wave 
equation 10 with boundary conditions 12 can be expressed as 
a finite sum of weighted normal modes. 


h h h (29) 

P(t, x,y,z)=YjYjYj 2* 1*2*3 b)</v 2 *3 * $ 

*1 k 2 =0 k 3 =0 


where ty kl k 2 k (x, y, z) are the normal modes of the homoge- 
neous part of equation 10 with boundary conditions (equation 
12). That is, 

y, z)=A klk 2 k 3 COS l k{ X COS 1*^ COS 1*3 2 (k 1; 

k 2 ,k 3 =0, 1,2, . . .) (30) 

where % k =k 2 ixJ\ y , ^ =k 3 it/l z and q* lMj (t) are the 

modal coordinates. Substituting equations 26 and 29 into 
wave equation 10, we get 


l l l 2 l 3 (31) 

Z Z Z [^*1*2*3 y < *»*1*2*3 (t) + 

4^=042=0^3=0 

4(^1 + fk 2 + 4 3 ) 2 ‘ ^* 1 * 2 * 3 ^ y, £)?4 1 4 2 4 3 W] = 

Pi P2 

4po£ Yj QrmWnmix - a, y - c)[H(x - a) - H{x - b)][H(y -c)- 

n=l m= 1 


H(y-dm Z ) 


Multiplying equation 31 by ^^^(x, y, z), integrating the 
equation over entire volume of the enclosure 400, using the 
orthogonality condition, and adding the modal damping term 
to the equation, we obtain 


#y4j4 2 43 0) "*■ 2^4 1 ^>k 2 %k 3 (L> kLk 2 k 3 ^kLk 2 k 3 (0 w 4j4 2 43 ?4j4 2 43 (0 


(32) 


"irZZ ^1*2*3 nmq^t) 


where are modal damping ratios, 


u kl k 2 k 3 = Coyj^ +%l 2 
V = ljyl z 


® kLk 2 k 3 nm 



^k Y k 2 k 3 {xL +a, yi +c, 0 )<p nm {xL, yi)dS { 


65 
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Now that the solution for acoustic pressure has been 
described, the state- space models of the system shall be 
described. Let the state variables be defined as 

*=(*>, y n m p ) T &w“, ( 33 ) 

where x Mm p =q Mm , y„ m p =q,, m , x JM ^ a ^q Wj , y* 1 * 2 * 3 “=q* ] * 2 * 3 , 
n=l, 2, . . . p : ; m=l, 2, . . . p 2 ; 1^=0, 1, . . . , l x ; k 2 =0, 1, . . . , 
1 2 ; and k 3 =0, 1 , . . . , 1 3 . If the external force can be expressed 
as f(t, x 1; y 1 )=f(t)g(x 1 , y 1 ), f wm (t) can be expressed as 



14 


-continued 

0 

0 

011 (*11, >ll) 

■■■ 01l(*lr, yir) 

0 

0 

<Pp lP2 (xn, J11) 

■■■ <f>p lP2 (xi r , yir) 




g(x 1, yO<Pnm(Xl,yi)dS l \f(t) 


= 7rmf(t) 


Then equations 28 and 32 can be re-written as 


yL = ~ M L x L - ^nm^nmyL + 7 ^ Y <Pnm(Xli, Vl/)«/(0 + T^/W 


P2 / 2 

"^1*2*3 — / I y ( Xk^k2^nm ynm ^1*2*3 — ~ ^1*2*3 _ ^4l^2^3 ^1*2^3 


In the state-space form, the model for a 3-D enclosure If we place k microphones in the position of (x., y ., z.), 

under consideration is given by i=l, . . . , k, the outputs of the system can be written as 

P(t)=(p(t, X, , y, , Zi p(t, x*, y zYi) T . Since matrix E is 
Ex(t)=Ax(t)+Bu(t)+Df[t) (36) . V } K *\ 9 19 J 19 u ’ 9 FV 5 k)) . 

35 invertible, the state-space model (equation 36) can be re- 

Matrices E, A, B, and D are given by written as 


[£ii 0 1 A n 0 

E = A = 

£21 £22 J [0 A 22 


Jc(r) = ET l Ax(t) + E~ l Bu(t ) + E~ l Df(t) 
Pit ) = Cx(t ) 


hpo Si [ 0 


where 

C = [C[C' 2 -C k ] T 


where, E n =I and A 11 =diag(A 11 ” m ) are square matrices of 

order p 1 p 2 , E 22 =I and A 22 =diag(A 22 A:i * 2 * 3 ) are square matrices The C matrix has k rows, and C. has the form 

of order (l 1 +l)(l 2 +l)(l 3 +l),B 11 isap 1 p 2 xr matrix, and D t t is 

a pip 2 x 1 matrix. Moreover, matrices E 21 , A 11; A 22 , B 11; and 50 * •’ =1 ’ \ 001 X ^ v Zz ’ ’ ’ 71/3/3 Zl ’ 

D ll are given as follows 


0 0 ■■■0 0 

0 arooin "'0 aooi Pl p 2 


Now that the state space model has been described, the 
transfer function model shall be described. For single-input 
single-output (SISO) systems, if we let f(t, x 1? yJ^O and take 
55 the Laplace transform of equations 28 and 32, we get 


0 0 ■■■0 0 

.0 a^n ■■■0 a W3PlP2 

0 1 
~OJ nm — 2^ nm OJ nm 


, , 1 (39) 

(S Z + 2 ZnmWmS + ^in)Qnm(^) = j^j^tnmiXu , Jll)^^) 


^1*2*3^ ^*1 *2*3 )^*1*2*3 (^ — 


Z Z “*1*2*3 w 


-Yk^k 3 - 24 ^ 2 * 3 ^ 1 *: 
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-continued 


Qk i k 2*3 CO : 


Js z 


(41) 


O ^£*1*2*3 W * 1*2*3 S 1*2*3 


X’' 


(r)V 2 P(r, s)dV - 


(47) 


^4- — *2*3wn$nm(*ll, Jll) 


ZjZj ^2 

n=l m=l 




hwL 


f/(5) 


5 2 P 

— %-(r)P(r, j)rfV 

c 0 Jv 


Pl P2 

~ -p 0 5 2 ^ 2 7>«£WO 

n=l m=l 


where 


where 

10 


4pq 

Vhp p S i ' 


^ = £%( r ) 

<Pnm(x -a, y- c)[H(x- a) - H(x - b)][H{y - c) - H(y - d)]6(z)dV. 


Taking the Laplace transform of equation 29, we can obtain a Using Green’ s first formula in terms of the appropriate vari- 

SISO transfer function for the output being the acoustic pres- a bles, the first term in the left hand side of equation 47 can he 

sure measured at the microphone located at (x 0 , y 0 , z 0 ) as 2 o expressed as 


. 


*1=0 *2=0 *3=0 


0*!*2* 3 (*O. JO, Zo) 


O ■*" ^£*1 *2 *3 W * 1*2*3 5 W * 1 * 2*3 


zz 


^kik2k^nm < Pnmi.^ll , Jll) 
S 2 + 2^ nm (O nm S + (0 2 m 


(42) 


x 


%(r)V 2 P(r, s)<jfV = 


(48) 


J " V j(r)V P(r, s) • «(r)af.s - V tyj(r) ■ V P(r, j)rf V 


We can then calculate 


To obtain the transfer function of the enclosure (i.e., equa- 
tion 42), we used “hard” boundary conditions (i.e., equation 
11 or equation 12. For more general boundary conditions, 
such as the soft boundaries, the following transfer function is 
used. On the boundary S, there exists the following relation- 
ship 


P{r, s)=Z(k s)U{r, s)-n(r ) (43) 

where r=(x, y, z), P(r, s) and U(s) are the Laplace transforms 
of p(t, r) and v(t, r) respectively, and Z(r, s) is the wall imped- 40 
ance. Substituting equation 43 into the Laplace transform of 
equation 9, we can get 


r 1 

%-(r)VP(r, s)-n(r)dS = - p 0 sj ' [Z^G-fr)] 
J s £i 

£v't J (r)-VP( r , s)dV = 2 fijiQiW 

£%{r)P(r, s)dV = Qj(s) 


where 


V P(r, s) ■ n(r) = -p 0 s 


P(r, s) 
Z(r, s) 


where r is S 


(44) 45 


Dji(s) 


■L 


Vj(sy9i(s) 

Z(r, s ) 


dS, p ri 


■£ 


V ¥•(/■) V %{r)dV. 


Assume that the solution of the Laplace transformed equa- 
tion of the wave equation (equation 1 0) can be written as 


Equation 47 becomes 


i 

P(r, s ) = Y i 'V,(rms) 

i = 1 


(45) 


55 


V [p 0 sDji(s) + pji]Qi(s) + S —Qj{s) 
£i c o 

y = l,2,- ,1 


pi p2 

= P0 S ^ jnm Qnm (0 

n= 1 m=l 


(49) 


where the following orthogonality condition are imposed on 
mode functions 


60 


Note that in practice it is often found that D^s) for j^i are 
negligible, in that case the expression for Q/s) can be written 
as 


x 


%(r)Vj(r)dV =6ij 


(46) 


Multiplying the Laplace transform of the wave equation 65 
equation 10 by 0 7 (r) and take the volume integral of the 
equation to get 


Qj(s) = j 2 — X2X £ Z 

— + posDjj(s) + fijj n=l m=l 


( 50 ) 
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Substituting equations 50 and 39 into equation 45, we get 
the transfer function for the enclosure with soft boundaries as 


V Po^4 %(r 0 ) 

Sb S L—i h Pp S i s 2 + pqCqsDh(s) + cl pa 

i=i 


' Pl P2 

zz 


VinmQnmix n, ^n) 
s 2 4- '2£ nm co nm s 4- co 2 f 


(51) 5 


10 


Now that the modeling has been described, the linearized 
dynamics of a structure-acoustic enclosure can be expressed 


in the state-space form as: 

15 

x(t)=Ax(t)+Bu(t)+B N N(t) 

(52) 

y(t)=Cx(i)+Du{i)+D N N{i) 

(53) 

z(t)=Cyx(t) 

(54) 20 


where x, u, y, and z are n-, m-, m-, and 1 -dimensional state 
vector, input vector, control-output vector, and performance- 
output vector, respectively (all functions of time t). The state 
vector x(t) consists of the acoustic state variables (pressures 2 5 
and rates) as well as structural vibration state variables 
(modal amplitudes and rates for the elastic modes). N(t) is the 
external disturbance vector consisting of external noise and/ 
or forces (e.g., due to the aircraft engine or air- friction at the 
airframe). The control -output y(t) is utilized for generating 30 
the feedback control signal u(t), while the performance output 
z(t) denotes a combination of noise levels (pressures) and 
vibration (accelerations) at locations where the performance 
is monitored. A, B, C, D, B^, are appropriately dimen- 
sioned matrices. The relationship between u(t) and y(t) can 35 
also be expressed as the transfer function 

P(s)=C(sI-A)- l B+D (55) 

The controller design is based on the concept of passivity. 

An mxm rational matrix G(s) is said to be positive real (PR) 40 
if 

(i) all elements of G(s) are analytic in Re(s)>0; 

(ii) G(s)+G*(s)>0 in Re(s)>0 (where G* denotes the con- 
jugate-transpose); or equivalently, 

(iia) poles on the imaginary axis are simple and have 45 
nonnegative-definite residues, and 

(iib) G(ja))+G*(jci))^0 for co€(-oo 5 00 ) 

For single-input, single-output (SISO) systems, condition 
(iib) is equivalent to: Re[G(j oo)]i^O, or equivalently, the phase 
of G(joo) remains between -90° and +90°. 50 

G{s) is said to be weakly strictly positive real (WSPR) if 
G(s) is stable (i.e., has all the Smith- McMillan poles in the 
open left-half of the complex plane), and the inequality in 
(iib) is strict (>), which requires the transmission zeros to be 
in the open left-half plane. The stability requirement can be 55 
weakened as follows to define the least restrictive class of 
positive real matrices, which permits imaginary -axis poles. 

An mxm rational matrix G(s) is said to be marginally 
strictly positive real (MSPR) if it is positive real, and 

G(/(o)+G*(/(d)> 0 for coe(-co, 00 ) 

Note that if G(s) is stable and MSPR, it is WSPR. 

The most important property of positive real systems is that 
they are robustly stabilized by any MSPR controller C(s) if 65 
none of the purely imaginary poles of G(s) is a transmission 
zero of C(s). Thus a PR system is robustly stabilized by any 


WSPR controller. The stability is guaranteed regardless of 
modeling inaccuracies and uncertainties. 

This robust stability property offers an attractive approach 
to robust controller synthesis for PR systems that can provide 
optimal performance along with guaranteed stability. How- 
ever, the transfer function of an acoustic -structure system is 
not positive-real, (i.e., it is a non-passive system). The 
approach of this invention is to “passify” (i.e., render passive) 
the system using compensation. Subsequently the feedback 
loop is closed with a suitably designed positive-real control- 
ler. 

The passifying compensators include series compensation, 
output feedback compensation, feedforward compensation, 
hybrid compensation (i.e., combination of series, feedback, 
and feedforward), and sensor blending and/or actuator allo- 
cation (applicable to systems having multiple sensors and 
actuators). FIGS. 4 a-e shows block diagrams of different 
types ofpassifications. FIG. 4 a illustrates series passification. 
FIG. 4b illustrates feedback passification. FIG. 4c illustrates 
feedforward passification. FIG. 4 d illustrates hybrid passifi- 
cation. FIG. 4e illustrates sensor/actuator blending. In each of 
FIGS. 4a-4e, P(s) indicates the non-passive plant and P c (s) 
indicates the passified plant. FIG. 4-e shows one of many 
possible configurations for hybrid passification. For single- 
input, single output (SISO) systems, passification can be 
accomplished in a fairly simple manner by using Bode plots, 
whereas for multi-input multi-output (MIMO) systems, pas- 
sification is more complex and is best performed in state- 
space setting using tools such as linear matrix inequalities. 

Different types of non-passive systems warrant different 
types or combination of different types of passification meth- 
ods. An important consideration in the design of a passifying 
compensator is the robustness of passification. The stability 
robustness depends on the robustness of passification as the 
passivity is not inherent to the system, but depends on the 
passifying compensator. The structural-acoustic system is 
usually non-minimum-phase. Therefore the approach of this 
invention is to first render the system minimum-phase by 
using a feedforward compensation Cff(s). For many struc- 
ture-acoustic systems, a constant feedforward matrix Dff is 
sufficient. The design of such a matrix can be accomplished 
by using linear-matrix-inequalities (LMI). Only the mini- 
mum necessary feedforward gain should be used for passifi- 
cation because it tends to limit the effective feedback gain and 
hence the performance. 

Designing a series passifying compensator can be accom- 
plished by using LMI techniques. Alternatively, if large num- 
bers of actuators and control sensors are available, the tech- 
niques of optimal blending of sensor signals and optimal 
allocation of control inputs are used. After designing the 
passifying compensator, the passivity of the compensated 
plant is verified using experimentally obtained frequency 
response data of the system and numerically generated fre- 
quency response data of the passifying compensator. In addi- 
tion, numerical computations are also performed to verify 
that the system transfer function remains ‘robustly PR’ in 
spite of parameter variations within a specified range (such as 
+/-25%). If the system does not remain ‘robustly PR’, the 
passification is redesigned to accomplish robust positive real- 
ness. 

Designing an optimal controller that is also PR. The meth- 
ods described in the paper “Synthesis of LQ-Optimal Con- 
stant-Gain Positive-Real Controllers”, Control and Intelli- 
gent Systems, 29(3): 65 -73, 2001 by A. G. Kelkar,Y. Mao, and 
S. M. Joshi, hereby incorporated by reference in its entirety, 
for optimal constant-gain positive real controllers are used. 
Closed-loop simulation studies are performed, using a white 
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noise disturbance input signal oo(t), to determine the reduc- 
tion in noise and vibration. If the reduction is not adequate, an 
optimal PR linear-quadratic-Gaussian (LQG) controller is 
designed. The performance is verified via simulation. If the 
performance i s not adequate, additional actuators and control 5 
sensors are installed and the design process is repeated. If the 
performance is satisfactory, the overall controller, consisting 
of the passifying compensator and the PR feedback control- 
ler, is realized as a digital filter and is implemented as the 
“energy-extracting controller” shown in FIG. 1. 10 

Now that the theory has been described, the details of the 
passification techniques shall be described. Passification can 
be performed via series compensation. Series compensation 
can be performed either at the input (pre-compensation), as 
shown in FIG. 4(a), or at the output (post-compensation). We 15 
shall only present the pre-compensation technique. The series 
post-compensation is similar. 

Consider a non-passive plant: 

20 

( x = Ax + Bu' (56) 

G(s)~< 

( y = Cx + Du' 


Series passification can be used only for Lyapunov- stable, 
minimum-phase systems. Acoustic- structure systems are 
always asymptotically stable, but they can be non-minimum- 
phase (i.e., one or more transmission zeros of the realization 
[A, B, C, D] can be in the closed right half of the complex 30 
plain. If the system is non-minimum-phase, it must first be 
rendered minimum phase via feedforward compensation. If 
an mxm constant feedforward gain matrix Dff is used, the 
resulting system is [A, B, C, (D+D^)] , and its inverse is given 
by: [(A-BD^C), BD e /\ D ejr -1 C, D ejr -1 ] where D^D+D^ 35 
is assumed to be nonsingular. Therefore the transmission 
zeros of [A, B, C, (D+D^)] are the eigenvalues of 
(A-BD ejr -1 C). Thus, if the system [A, B, C, D] can be stabi- 
lized by constant-gain output feedback F, it can be rendered 
minimum-phase by a constant feedforward gain gain 40 
which is obtained from: F=(D+D^) _1 . F can be obtained by 
using multivariable pole placement methods or by trial and 
error by setting Dj to be a positive diagonal matrix and 
increasing its elements until the inverse system is stable. It is 
advisable to find the smallest such that the inverse system 45 
is stable. The resulting minimum-phase system is henceforth 
denoted as: [A, B, C, D]. 

The series passifier (FIG. 4(a)) is given by: 

50 

( x c = A c x c + B c u (57) 

C 5 (s)~I 

\u'=C c x c + D c u 


Then the compensated system is given by: 


For the compensated system (G c (s)) to be passive, we need 
(A, B, C, D) to satisfy the following inequality 


a t p+pa 

PB 

+ 

c d 

r ro -r 

c b 

b t p 

0 


.0 / _ 

[-/ 0 

_ o / _ 


Matrix P is a symmetric positive definite solution to the 
inequality (i.e., Equation 60). Existence of such P ensures 
positive realness of Gc(s). Note that Equation 60 is nonlinear 
in matrix variables A c , B c , C c , D c , and P, and therefore, does 
not constitute an LMI. This nonlinear matrix inequality, how- 
ever, can be converted into an equivalent LMI using a judi- 
ciously chosen nonlinear matrix transformation. Given below 
is the procedure to convert this nonlinear matrix inequality 
into an LMI. 

Let us suppose matrices P and P -1 can be partitioned as: 


P = 


Y N 

n t * 


r l 


x M 
m t * 


(61) 


where X and Y are symmetric and full rank matrices with the 
restriction that 

XY+MN r =I (62) 

Further, if we define 


rii 



/ Y 

0 n t 


(63) 


then, since PP 1= I, ll 1 and II 2 have the following properties 


prii = IE and nfprii = n[n 2 


x 

i 


i 

Y 


(64) 


Let 


n = 


o 

/ 


[ x = Ax + Bu 
[ y = Cx + Du 


Now pre- and post-multiplying Equation 60 by U T and IT, 
(58) respectively, does not change the inequality, i.e., 

60 


where 




nf^p+p^ri! n \pb 


o -nic r 

r * i r a bc c 1 _ \bd c ] 

(59) 

+ 

-crii -(d + d t ) 

U’Ho J’*=LJ 

(p r p)ni o 



( 65 ) 


C[C DC c ], and D = DD C . 


65 


where 
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-continued 

n[ pa rii = nlAYh 


0 

A 

BC c 

X 

I 

N _ 

.0 

Ac . 

m t 

0 


ax + bc c m t a 

YAX + YBC c M t + NA c M t YA 


Now consider the following change of variables 

A = YAX+ YB C c M r +NA c M t 

B=YBD C +NB C 

C=C c M r 

D=D C 

Using transformation equations (66), we get 


nfPArii = 


AX +BC A 
A t KA 


(AX + BC) A t 


n\PB = niB = 


cn { = [C DC c ] 






, \‘ 

°1 

\ BD A 

BD 


n\ 

u \- 

B 


X I 
M t 0 


= \CX + DC C\ 
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Passification via feedforward (parallel) compensation shall 
now be described. Consider a plant: 


5 



x = Ax + Bu' 
y = Cx + Du ' 


(69) 


and feedforward compensator (FIG. 4 (c)): 


15 




{ 


x c = A c x c + B c u 
yi = C c x c + D c u 


(70) 


Then the compensated system is given by: 


( 66 ) 


20 


(67) 



f x = Ax + Bu 

G e (s)~\ 

\y -Cx + Du 

where 


\x 

1 ~ [/i o 

x=\ 

L A = 1 

U 

J [0 A c 


,B = 


and 


C=[C C c ], D = D + D c . 


(71) 


(72) 


30 For the compensated system to be passive, we need, as in the 
previous cases, 


a t p+pa 

PB 

+ 

c 

b 

T \° 

-i 

c 

b 

b t p 

0 


. 0 

/ _ 

[-/ 

0 

_ 0 

/ 


Substituting Equations 67 in Equation 65, we get an LMI in 
variables X, Y, A, B, C, and D. 

Pre- and post-multiplying equation 72 by IT 7 and IT, respec- 
40 tively, we get 


A** C) O 

A + A t YA + A t Y C) 

d t b t -cx-dc b t -c d** 


(68) 


45 


where A* *=AX+XA r +BC+C r B r , D ,| " |, =-DD-D 7 b 7 ’and (*) 
indicates the transposes of their corresponding off-diagonal 
counterparts. 

If the solution to Equation 68 is feasible, the passifying 50 
compensator Cs(s) can be recovered by following the step- 
by-step procedure below. 

1) Solve equation 68 to obtain X, Y, A, B, C, and D 

2) Construct M, N, and P such that they satisfy Equation 64 . 5 5 
(Note that M and N should satisfy Equation 62. Also, 
since 


n[(A r p+PA)ni n \pb 
(Fpjri! o 



-nic 7 
-(d + d t ) 


<o 


where 


nfpAn: =nlAYh = 


AX A 

YAX + NA c M t YA 


Now consider the following change of variables 

a=yax+na c m t 


B=NB C 


C=C c M r 


(73) 


(74) 


X I 1 

>0, Y >0, 

/ Y 


60 


D=D, 


Using equations 75, we get 


(75) 


and X-Y ^0, 1-XY is nonsingular, i.e., we can always find r ax a 

nonsingular square matrices M and N satisfying equation 62) 65 nfPArii = , ^ 

3) Solve equations 66 in reverse order to obtain Ac, Be, Cc, 
and Dc. 


( 76 ) 
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-continued 


n[Y PFii = 


n [pb= n‘ 2 B = \ 


cn { = lex + c c 


p z - and p z , i=l , 2, . . . k, denote the lower and upper bounds on 
p z .. Define 


Pi = ^(Pi +£ f .) and (5; = / = 1, 2, ...* 


Expanding p z - as p z =p z +a z 8 z where a z €[-l, 1], the system 
matrices can be written as: 


Substituting equations 76 in equation 73 we get an LMI in 
variables X, Y, A, B, C, and D. 


A(p) = A(jn) + y aiSiAr, B(p) = B(p) + £ a-.S-.B, 


AX + XA 1 (*) (*) 

A + A T YA+A t Y (*) <0 


\B‘ -CX -C B'Y + B-C D" 


where D _L =-(D+D r +D+D 7 ) and (*) represents elements 
which are transpose of their corresponding off-diagonal 
counterparts. The design process will again involve the steps 
given in previous paragraphs. In this case also, similar to 
series passification, it is required that the non-passive plant be 
open-loop stable. 

Constant-gain feedforward passification shall now be 
described. The simplest feedforward compensation is a con- 
stant gain Dc rather than a dynamic compensator. In that case, 
the compensated system is PR if and only if there exist matri- 
ces P=P r >0 and Dc such that 


B‘P-C -{D + d c ) -(D + D c y 


C(p) = C(p) + y a,S,C,; D(p) = D(p) + y <*,<$, D, 


where n=[ji 1; 1 1 2 , • • • H*] r and A„ B„ C„ D, are constant 
matrices. The ith element (i=l, 2, ... k) of the jth vertex p 7 , 
(j=l, 2, ... 2*) of (3 is given by: 

phPi+YiA (Y ij=- 1 or 1 ) ( 84 ) 

Then the system matrix at p^p 7 is given by 


A(p J ) = A{p) + JijdjAi 


B^), C(p 7 ), and D(p 7 ) are similarly defined. 

The following theorem presents a sufficient condition for 
the system to be robustly passive V^€B. Note that passivity is 
defined only for ‘square * systems (i.e., with the same number 
of inputs and outputs). Suppose l=r. The system S(p) is pas- 
sive if there exists a symmetric positive-definite matrix P that 
is a solution of the following system of 2 k LMIs: 


Dc is constrained to be symmetric and non-negative definite 
and a solution that minimizes the trace of Dc is obtained. 

Passification via sensor blending and control allocation 
shall now be described. In some cases, many sensors (accel- 45 
erometers and/or microphones) may be available for measur- 
ing the vibration accelerations and pressures at various loca- 
tions. In such cases, sensor blending can be performed 
(matrix M in FIG. 4(e)) to passity the system. Similarly, if 
many actuators (piezos and/or loud speakers) are available, 50 
control allocation can be performed (matrix N in FIG. 4(e)). 
These passification techniques for obtaining a system that is 
robustly positive real in the presence of parameter uncertain- 
ties are described herein. 

55 

Consider the system S (p) given by: 

x=A(p)x+B(p)u; y=C(p)x+D(p)u (79) 

where peR k denotes the vector of uncertain parameters and x, 
u, y are n-, r-, and 1 -dimensional state, input, and output 60 
vectors respectively. A(p), B(p), C(p), D(p) are appropriately 
dimensioned matrices that are assumed to be affine functions 
of the parameter vector p. The parameters p z (i=l, 2, . . . k), 
which are components of p, are assumed to lie in a hyper- 
rectangular box in the parameter space, defined as: 


" A T (p J )P + PA(p j ) PB(p J ) - C T (p J ) 

_ B T (pj)P-C(pj ) ~(D(pj) + D T (pj)) _ 


for y=l,2, 


For sensor blending where l>r, it is required to find a 
blending matrix MeR rx/ such that the system S^=[A, B, MC, 
MD] is passive. Further, it is required to maximize the size of 
the parameter box (3 in which the system remains passive. 
Consider the parameter box defined by {p :p z € [u z -0 6 Z , u z +0 § z ] , 
i=l, 2, . . . k} where 0 is a dilation parameter. The ith element 
of the vertex p 7 , j : 1,2,... 2 k , of the parameter box with 
dilation 0 is given by: 

pl=\Wfih i (i=l, 2,...k) (86) 

The problem is to find a blending matrix M which will maxi- 
mize the dilation 0 while passifying the system. 

The system matrices at the vertices have the form: 


A{p i ) = AQi) + VjSiAi ; B(pj) = B{p) + 9^ vABi 


P={P : P, e [n^PT i= l,2,...k} 


(80) 
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-continued 

k k 

C(pj) = C(p) + r,ACi\ D(pj ) = D(p) + (>£ y^D, 

i = 1 i'=l 


( 88 ) 

5 



Substituting equations 87-91 in 85, the optimal sensor 
blending problem can be stated as follows: Find matrices 
PeR”*”, (P=P r >0) and MeR r * z which will maximize 0 subject 
to: 25 


A T sj P+PA sj PB sj -C T sj M T 1 (92) 

B t s .P-MC sj -(MD sj + D t s .M t 

30 

forj=l,2, ...2*. 

After sensor blending, the resulting pair (MC, A) must 
retain observability. The following constraints are added to 35 
ensure observability for p=g: 


ZQu, P, M) < 0 


' A r (ji)P + PA(ji) PB(fi) - C T (ft)M T (95) 

B t {ji)P- MCiji) ~(D + D t ) 

\A T sj P+PA sj PB sj -C T sj M T 
[ BljP-MCsj 0 

for j=l, 2, . . . 2 k . 

If the D matrix is added to a strictly proper system and K(s) 
denotes the feedback controller for the passified system (A, 
B, MC, D), the effective controller (from y=Cx to u) becomes 
K eif (s)=[I+K(s)D] -1 K(s)M. If K is a constant positive definite 
matrix, K eJf — >D _1 M as K is increased. Thus D has the effect of 
reducing and limiting the effective control gain. Therefore, if 
D is a design variable, it is desirable to make it as small as 
possible, which can be accomplished by placing an upper 
bound on it (i.e., D+D r <d m J). 

Note that even for the non- strictly -proper case, an addi- 
tional variable D can be used for making (A, B, MC, MD+D) 
passive. This usually increases the region of robust passivity 
but reduces the effective controller gain. 

For optimal control allocation (r>l), it is required to find a 
control allocation matrix NeR^ such that the system [A, 
BN, C, DN] is passive. It is also required to maximize the size 
of the parameter box (3 within which the system remains 
passive. Define 

[ A(pt)L + IA T Qi) B(n)N -ZC T (/u) 1 (96) 

Z(j*,Z,N):= 

B t (h) - C{p)Z ~(D(p)N + N T D T (fi )) 


~A t ( M )Q 0 +Q 0 A( M ) 

< 0 

(93) 

MC(p) 

-I 


Qo >£o 


(94) 


The problem can be stated as follows: Find matrices 
2eR" x ", (2=2 r >0) and NeR 7 ^ which will maximize 0 subject 
to: 


where € 0 is the specified lower bound on Q 0 , a measure of 
observability. 

Equations 92-94 define a generalized eigenvalue problem 
which can be solved by using LMI methods. If p z and p z are not 
known apriori, but only their nominal values (p z ) are known, 
it would be necessary to guess the corresponding 5 Z . and the 
maximum 0 obtained will depend on the choice of 5 z ’s. 

Note that, if D=0, the off-diagonal blocks in the passivity 
LMI (equation 85) must be zero. If uncertain parameters are 
present only in the matrix, then robust passification is pos- 
sible even when D=0. If, however, uncertain parameters are 
present in B and/or C, it will be necessary to have a nonzero 
D matrix in order to satisfy the sufficient conditions of equa- 
tion 92. Thus if the system to be passified by sensor blending 
is strictly proper, an additional unknown variable D will have 
to be included in order to make the system (A, B, MC, D) 
robustly passive. The generalized eigenvalue problem is re- 
stated as follows for the strictly proper case. 

For the strictly proper case, the solution is to find matrices 
PeR"«", (P=P 7 ’>0), DeR^, and MeR^ which will maximize 0 
subject to: 


„ A sJ L + ZA t s j B S jN - ZClj (97) 

Z(p,Z,N)<9 . . 

[N T Blj-C s jZ -(D sJ N + N T Dl)\ 

for j=l, 2, ... 2* 

A lower bound constraint on a measure of controllability 
Q c is also needed: 

T A(fi)Q c + Q c A T (fx) B(m)N (98) 

r r <0 

N T B T (fi ) -I 

Qc>s c l (99) 

where € c is the specified lower bound on Qc. 

When redundant actuators and sensors are present, it would 
be advantageous to perform both control allocation and sen- 
sor blending. If we let m^min (r, 1). The problem is to find 
matrices MeR mx/ , NeR rx ™ such that the system (A, BN, MC, 
MDN) is robustly passive (i.e., find a positive definite matrix 
P and matrices M, N such that 0 is maximized subject to 
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' A r {fi)P+PA{fi) PB(pi)N -C T (fi)M T 1 (100) 

N T B T (&P - MC(v) ~(MD(jj.)N + N T D T (fX)M T 

[ A T sj P + PA sj PB sJ N - C T sj M T 
9 N T B T sj P - MC s j -(MD s jN + N T D T sj M T 

for j=l, 2, . . . 2* or, using the dual LMI for passivity, the 
problem is to find a positive definite matrix and 2 and matri- 
ces M, N such that 0 is maximized subject to 

' A{pL)Z + ?A T Ui) B(fi)N 1 (101) 

N T B T {ji) - MC(jd)Z +N t D t (h)M t ) 

+ A T sj 2 B sj N - SCj-M 7 

^ N T B T sj -MC sj 2 -(MD sj N + N T DljM T 

for j=l, 2, . . . 2*. Unfortunately, quadratic terms involving 
products of (M, N), (P, N) and (2, M) are present in equations 
100 and 101, and are therefore no longer LMIs. However, a 
possible procedure for solving this problem is to iterate 
between sensor blending and control allocation. After select- 
ing an m ^ min (r, 1), a full rank M matrix is arbitrarily selected 
and the optimal control allocation problem is solved, yielding 
an N matrix. Next, the optimal sensor blending problem is 
solved for this fixed N, yielding a new M. The procedure is 
stopped when the 0s obtained are sufficiently close and large. 

The techniques of passification have been described. FIG. 
5 shows a passivity based controller structure 400 for a one 
dimensional enclosure system 402 using the techniques 
described above. Speaker 404 is the control speaker and 
speaker 406 is the disturbance-generating speaker. C^(s) 
compensator 408 denotes the transfer function of the feedfor- 
ward part of the controller 400. C 5 (s) compensator 410 
denotes the transfer function of the series part of the controller 
400. C^(s) compensator 412 denotes the transfer function of 
the feedback part of the controller 400. N(s) and y per /s) 
denote the Laplace transforms of the disturbance noise and 
the sensed noise respectively. 

FIG. 6 shows the experimentally obtained frequency-do- 
main noise suppression response for the enclosure 402 for a 
broadband disturbance noise. FIG. 7 shows noise suppression 
in a three-dimensional reverberant structure for a three-tone 
disturbance noise at frequencies of 100 Hz, 225 Hz, and 290 
Hz that was obtained via simulation. FIG. 8 shows the experi- 
mentally obtained noise suppression response for the struc- 
ture of FIG. 7. The suppressed noise levels correspond to the 
“controller on” condition. 

Now that the modeling, passification, and controller design 
has been described, an example shall be provided. Analyti- 
cally derived models (i.e., equations 37 and 42), for the acous- 
tic-structure interaction dynamics in 3-D enclosures having 
configuration shown in FIG. 1, are used to obtain the finite- 
dimensional models for the laboratory apparatus shown in 
FIG. 9. The laboratory apparatus 500 is a rectangular wooden 
box with one of its surfaces being the aluminum plate. This 
surface can be set into vibrations to generate an acoustic 
disturbance inside the box using a piezoactuator attached at 
the center of the aluminum plate. System parameters for this 
3-D acoustic enclosure are given in Table 1. 


TABLE 1 


Parameter 

Value 

Parameter 

Value 

l x 

0.3048 m 

V 

0.3810 m 

\ z 

0.7874 m 

a 

0.0063 m 

b 

0.2985 m 

c 

0.0063 m 

d 

0.3747 m 

h 

0.000813 m 

E 

71 * 10 9 N/m 2 

P 

0.33 

P 

2810 kg/m 3 

c 0 

343 m/s 

Po 

1.13 kg/m 3 

x o 

0.1270 m 

yo 

0.1524 m 

zo 

0.3112 m 

x i 

0.1524 m 

Yi 

0.1905 m 


Using the parameter values of Table 1 and assumed modes 
approximation, a transfer function model of the system is 
obtained using equation 42. For finite dimensional approxi- 
mation, the first five modes in each x, y, and z direction are 
considered, which yield a total of 125 modes. 

In order to validate the analytical model, an experimental 
frequency response was obtained using swept sine excitation 
of the piezoactuator as well as the speaker and measuring the 
acoustic pressure at the sensing microphone. The sensing 
microphone is located at location (171, 177, 285) mm. Using 
this experimental frequency response data, the state-space 
model of the system was derived using system identification 
techniques. The system identification toolbox SOCIT (devel- 
oped at NASA LARC) was used for this purpose. The SOCIT 
toolbox uses ERA (Eigensystem Realization Algorithm) 
method, which is based on Markov parameters and singular 
value decomposition techniques. The ERA uses modal ampli- 
tude coherence for ranking of the most effective modes of the 
system. For the system under consideration, this identifica- 
tion algorithm gave a very satisfactory match. The identified 
model is a 40th order model. It is to be noted that identifica- 
tion becomes increasingly difficult as we go from SI SO to 
MIMO systems. For the MIMO case, it is very hard to obtain 
a good match and it invariably yields very high order modes. 

The natural frequencies obtained from the identified sys- 
tem model were compared with the natural frequencies deter- 
mined from analytical model to assess the accuracy of the 
analytical model. The comparison is given in Table 2. The 
symbol * in the second column denotes the frequencies that 
are missing from the analytical model but appear in identified 
model. Similarly, the symbol # in the third column denotes the 
frequencies that are missing in the identified model but appear 
in analytical model. The missing dynamics in the analytical 
model is attributed to the fact that the analytical model does 
not account for actuator, sensor, and other hardware dynamics 
and the boundary conditions used to derive analytical model 
are only approximate. In particular, the modes 2, 7, 11, 14, 
and 19 are missing in the analytical model. The modes 3, 9, 
10, 18, 21, 22, and 24 are missing in the identified model. The 
modes 9 and 10 are very close to 8 and with appropriate 
damping may collapse onto one. This may be the case for the 
missing modes 9 and 10 in the identified model. For the same 
reason, modes 18 and 24 are missing in the identified model. 
The modes 3, 21 and 22 seem to be highly damped in the 
experimental data which is why the identification algorithm 
cannot detect. 

The frequencies listed in Table 2 do not account for bound- 
ary impedance, i.e., it assumes the perfectly hard boundary 
with zero transmissibility. For controller design, the identi- 
fied system model is used as the design model since it repre- 
sents the system more accurately than the analytical model. 
Although the natural frequencies from analytical model 
match reasonably well with natural frequencies from experi- 
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mental data, the damping on each mode needs to be adjusted 
manually since there is no systematic method to derive it 


TABLE 2 



Frequency (Hz) 

Mode no. 

Analytical 

Experimental 

1 

37.5 

39.2 

2 

* 

64.0 

3 

80.8 

# 

4 

106.4 

98.4 

5 

149.8 

125.5 

6 

153.2 

157.7 

7 

* 

182.4 

8 

217.8 

222.4 

9 

221.4 

# 

10 

222.1 

# 

11 

* 

239.0 

12 

254.4 

251.3 

13 

264.8 

262.1 

14 

* 

280.5 

15 

323.4 

327.1 

16 

337.1 

351.0 

17 

382.3 

381.0 

18 

384.6 

# 

19 

* 

416.7 

20 

425.7 

423.5 

21 

435.6 

# 

22 

438.4 

# 

23 

450.1 

454.5 

24 

453.6 

# 

25 

498.0 

479.0 


AAs can be seen from the above description, the approach 
taken here is to first attempt designing a feedforward com- 
pensation to passify the system. If this is not acceptable (for 
example, in the case of non-minimum-phase systems, large 
gain matrices may be required), a constant-gain feedforward 35 
compensation is designed to only to render the system mini- 
mum-phase. The resulting system is then rendered PR by a 
series (pre- or post-) compensation, or if applicable, by sensor 
blending and/or control allocation. The next step is to design 
a PR controller so that the resulting closed-loop response 40 
shows the desired noise reduction. 

All references, including publications, patent applications, 
and patents, cited herein are hereby incorporated by reference 
to the same extent as if each reference were individually and 45 
specifically indicated to be incorporated by reference and 
were set forth in its entirety herein. 

The use of the terms “a” and “an” and “the” and similar 
referents in the context of describing the invention (especially 
in the context of the following claims) are to be construed to 50 
cover both the singular and the plural, unless otherwise indi- 
cated herein or clearly contradicted by context. The terms 
“comprising,” “having,” “including,” and “containing” are to 
be construed as open-ended terms (i.e., meaning “including, 
but not limited to”) unless otherwise noted. Recitation of 55 
ranges of values herein are merely intended to serve as a 
shorthand method of referring individually to each separate 
value falling within the range, unless otherwise indicated 
herein, and each separate value is incorporated into the speci- 
fication as if it were individually recited herein. All methods 60 
described herein can be performed in any suitable order 
unless otherwise indicated herein or otherwise clearly con- 
tradicted by context. The use of any and all examples, or 
exemplary language (e.g., “such as”) provided herein, is 
intended merely to better illuminate the invention and does 65 
not pose a limitation on the scope of the invention unless 
otherwise claimed. No language in the specification should be 


30 

construed as indicating any non-claimed element as essential 
to the practice of the invention. 

Preferred embodiments of this invention are described 
herein, including the best mode known to the inventors for 
carrying out the invention. Variations of those preferred 
embodiments may become apparent to those of ordinary skill 
in the art upon reading the foregoing description. The inven- 
tors expect skilled artisans to employ such variations as 
appropriate, and the inventors intend for the invention to be 
practiced otherwise than as specifically described herein. 
Accordingly, this invention includes all modifications and 
equivalents of the subject matter recited in the claims 
appended hereto as permitted by applicable law. Moreover, 
any combination of the above-described elements in all pos- 
sible variations thereof is encompassed by the invention 
unless otherwise indicated herein or otherwise clearly con- 
tradicted by context. 

What is claimed is: 

1. A method for extracting acoustic energy and structural 
energy in an acoustic enclosure using a plurality of actuators 
comprising the steps of: 

obtaining a continuous -time multi-input multi-output 
state-space mathematical model of the acoustic enclo- 
sure; 

designing compensation to render the mathematical model 
passive in accordance with mathematical system theory 
if the mathematical model is not passive, thereby form- 
ing a compensated system that is passive; 

checking passivity of the compensated system; 

designing a passivity -based controller that extracts the 
acoustic energy and the structural energy such that a 
resulting closed-loop response provides a desired noise 
reduction; and 

extracting acoustic and structural energy from the acoustic 
enclosure using the actuators by controlling the actua- 
tors with the passivity-based controller. 

2 . The method of claim 1 wherein the step of obtaining a 
continuous-time multi-input multi-output state-space math- 
ematical model of the acoustic enclosure comprises the step 
of obtaining a mathematical model having the form according 
to the equation 

Ex(t)=Ax(t) +Bu(t) +D/[t) 

where A, B, D, and E are matrices given by 
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where E n =I and A 11 =diag(A 11 ” m ) are square matrices of 
order p 1 p 2 , E 22 =I and A 22 =diag(A 22 Ari * 2Ar3 ) are square matrices 
of order (1 1 + 1)(1 2 +1)(1 3 +1), B n is a p : p 2 x r matrix, D n is a 
p x p 2 x 1 matrix where matrices E 2 1 , A : : , A 22 ,B 1X , and D y 1 are 
given by 
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based controller includes designing a controller having a 
transfer function G^(s) wherein 


5 G sb (s) = 

Z POS 2 cl ^/(Po) V V ^ll) 
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where } ¥ i denotes the eigenmode function for the acoustic 
pressure expression obtained using the assumed modes 
method, r\ inm is the volume integral term consisting of inte- 
1 5 grand which is product of structural-acoustic eigenfunctions, 
Xxj is the damping of the ij -th structural mode shape, p 0 is fluid 
density at equilibrium, c 0 is the sound speed, S x is a boundary 
surface of the structure, h is a thickness of the enclosure, is 

the density of the plate, is the ij-th mode shape of the 
20 enclosure structure, and 


where h is a thickness of the enclosure, p 0 is fluid density at 
equilibrium, is a boundary surface of the structure, c 0 is the 
sound speed, V is the volume of the enclosure, ex’s are cou- 
pling coefficients describing the modal interaction between 
structural and acoustic modes, oo^ denotes natural frequency 
related to ij-th mode for the structure, ($ iJk denotes the acous- 
tical modal frequency for the ijk-th acoustic mode of the 
enclosure, . is the damping of the ij-th structural mode 
shape, t, iJk is the damping of the ijk-th acoustical mode shape, 
is the ij-th mode shape of the enclosure structure, and 
matrix D n indicate non-zero coefficients for the direct trans- 
mission terms which are functions of modal parameters. 

3. The method of claim 1 wherein the step of designing a 
passivity-based controller includes designing a controller 
having a transfer function G(s) wherein 
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where Z is the impedance. 
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5. The method of claim 1 wherein the step of designing 
compensation includes the step of designing a series passifier 
Cy(s) according to 
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40 wherein Ac, Be, Cc, and Dc are determined according to the 
steps comprising: 
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solving the equation 

A+A t 
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to obtain X, Y, A, B, C, and D\ 
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h is a thickness of the enclosure, p 0 is fluid density at equi- 
librium, Si is a boundary surface of the structure, c 0 is the 
sound speed, p^ is the density of the plate, (x, y, z) are 

normal modes of a non-homogeneous wave equation, 6Q 
with 1*,> 5/a, and 5 m being modal 
coordinates, t, iJk is the damping of the ijk-th acoustical mode 
shape, a’s are coupling coefficients describing the modal 
interaction between structural and acoustic modes, and is 
the damping of the ij-th structural mode shape. 65 

4. The method of claim 1 wherein the acoustic enclosure 
has a soft boundary and the step of designing a passivity- 


constructing matrices M, N, and P such that 
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and 

solving the equations A=YAX+YBC c M r +NA c M :r , 
B=YBD C +NB C , C=C c M r , and D=D C in reverse order to 
obtain Ac, Be, Cc, and Dc. 

6. The method of claim 1 wherein the step of designing 5 
compensation comprises the step of designing a feedforward 
compensator C^s) according to 
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wherein Ac, Be, Cc, and Dc are determined according to the 
steps comprising: 
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where D ± =-(D+D r +D+D+D r ) to obtain X,Y, A, B, C, and D; 
constructing matrices M, N, and P such that 
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Fill = n 2 and rh/tlli = 
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and 


solving the equations A=YAX+NA c M r , B=NB C , C=CJVl r , 
and D=D C in reverse order to obtain Ac, Be, Cc, and Dc. 35 

7. The method of claim 1 wherein the step of designing 
compensation to render the mathematical model passive com- 
prises the steps of: 

determining if a feedforward compensation will passify the 
system; 40 

if a feedforward compensation will not passify the system: 
designing a constant gain feedforw ard compensation to 
render the compensated system minimum-phase; and 
rendering the compensated system positive-real by at 
least one of series compensation, sensor-blending and 45 
control allocation. 

8. The method of claim 7 wherein the step of designing a 
passivity-based controller comprises the step of designing 
one of a dissipative linear-quadratic -Gaussian (LQG) type 
positive-real controller and a dissipative constant gain posi- 50 
tive-real controller. 

9. The method of claim 7 wherein the step of rendering the 
compensated system positive-real by at least one of series 
compensation, sensor-blending and control allocation com- 
prises the step of rendering the compensated system positive- 55 
real by at least one of series compensation, feedback com- 
pensation, hybrid compensation, and sensor-blending and 
control allocation. 

10. The method of claim 1 further comprising the step of 
redesigning the compensation if the passivity is not pre- 60 
served. 

11. The method of claim 1 further comprising the step of 
performing numerical simulations of the controller in the 
presence of a simulated broadband disturbance input. 

12. The method of claim 11 further comprising the step of 65 
redesigning the controller if the closed-loop response is not 
satisfactory. 
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13. A method for extracting acoustic energy and structural 
energy in an acoustic enclosure using a plurality of actuators 
comprising the steps of: 

obtaining a continuous -time multi -input multi-output 
state-space mathematical model of the acoustic enclo- 
sure; 

designing compensation to render the mathematical model 
passive in accordance with mathematical system theory 
if the mathematical model is not passive, thereby form- 
ing a compensated system that is passive; 
checking passivity of the compensated system; 
designing a passivity -based controller that extracts at least 
one of acoustic energy or structural energy such that a 
resulting closed-loop response provides a desired noise 
reduction; 

wherein the step of designing compensation comprises the 
step of performing sensor blending if there are redun- 
dant sensors; and 

extracting acoustic and structural energy from the acoustic 
enclosure using the actuators by controlling the actua- 
tors with the passivity-based controller. 

14. A method for extracting acoustic energy and structural 
energy in an acoustic enclosure using a plurality of actuators 
comprising the steps of: 

obtaining a continuous -time multi -input multi-output 
state-space mathematical model of the acoustic enclo- 
sure; 

designing compensation to render the mathematical model 
passive in accordance with mathematical system theory 
if the mathematical model is not passive, thereby form- 
ing a compensated system that is passive; 
checking passivity of the compensated system; 
designing a passivity -based controller that extracts at least 
one of acoustic energy or structural energy such that a 
resulting closed-loop response provides a desired noise 
reduction; 

wherein the step of designing compensation comprises the 
step of performing control allocation if there are redun- 
dant actuators; and 

extracting acoustic and structural energy from the acoustic 
enclosure using the actuators by controlling the actua- 
tors with the passivity-based controller. 

15. A method for extracting acoustic energy and structural 
energy in an acoustic enclosure using a plurality of actuators 
comprising the steps of: 

obtaining a continuous -time multi-input multi-output 
state-space mathematical model of the acoustic enclo- 
sure; 

designing compensation to render the mathematical model 
passive in accordance with mathematical system theory 
if the mathematical model is not passive, thereby form- 
ing a compensated system that is passive; 
checking passivity of the compensated system; 
designing a passivity -based controller that extracts at least 
one of acoustic energy or structural energy such that a 
resulting closed-loop response provides a desired noise 
reduction; 

wherein the step of designing compensation comprises the 
steps of: 

designing a constant gain feedforward compensation to 
render the compensated system minimum phase; 
rendering the compensated system positive-real by one of 
sensor-blending and control allocation; and 
extracting acoustic and structural energy from the acoustic 
enclosure using the actuators by controlling the actua- 
tors with the passivity-based controller. 



